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Abstract 

Akodon cursor occurs in dense rainforest from northern (8° S) to southern (26° S) states along the Atlantic coast of 
Brazil. Previous karyological and molecular data revealed two major clades, one including northern (8-1 5° S) and the 
other southern (19-26° S) populations. The center of geographic distribution (15-20° S), which included the state of 
Espirito Santo, was identified as a potential vicariance region. Since river barriers are among the most discussed 
models in the study of Neotropical diversification, we examined whether the Rio Doce (1 9° S) plays an important role 
in shaping the population genetic divergence of A. cursorby including samples from Espirito Santo in the analysis. 
Our results showed that the northern-southern division region in Atlantic forest was no coincidence with the presence 
of the Rio Doce by refuting the hypothesis that this river is an effective barrier to gene flow between populations. In- 
stead, we found evidence that isolation by geographical distance shaped the phylogeographical structure in the 
southern lineage. However, there is uncertainty about effectiveness of the processes involved and further studies 
based on wider sampling are needed. 
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Originally extending for 1,300,000 km 2 along the 
Brazilian coast and reaching into Paraguay and Argentina, 
the Atlantic forest has been reduced to less than 8% of its 
range (Gusmao Camara, 2003). This habitat loss has been 
closely linked to fragmentation but the region still harbors 
one of the highest percentages of endemic species in the 
world, with many species and even genera of vertebrates 
still being described (Rodrigues, 2005). Studies of the com- 
plex evolutionary history of this habitat have identified 
northern and southern components in the Atlantic Forest, 
with genetically structured populations and the differentia- 
tion of distinct groups throughout the Atlantic biome (Mus- 
trangi and Patton, 1997; Costa et al, 2000; Geise et al, 
200 1 ; Pellegrino et al. , 2005 ; Cabanne et al. , 2007; Noguei- 
ra and Fagundes, 2008). These findings suggest that a com- 
mon mechanism may have played a role in shaping the 
distributions of multiple taxa in this biome. Among several 
hypotheses for the diversification of rainforest biotas (Mo- 
ritz et al, 2000), evolution in palaeorefuges (Carnaval and 
Moritz, 2008; Thome et al. ,2010) and the influence of geo- 
graphical barriers (Pellegrino et al., 2005; Cabanne et al., 
2007) are the two most commonly invoked mechanisms. 
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The rodent Akodon cursor Winge, 1887 (Cricetidae) 
is widely distributed in rainforest along the Atlantic coast 
of Brazil (Musser and Carleton, 2005), from the northern 
states of Paraiba and Pernambuco (8° S) through Bahia, 
Espirito Santo, Rio de Janeiro, Minas Gerais and Sao Paulo 
to the northern region of Parana state in the south (26° S). 
This species has a peculiarly high karyotypic variability, 
with the diploid number varying from 2n = 14 to 16 and 
fundamental numbers (FN) ranging from 18 to 26 (Fagun- 
des et al, 1998). Although there is no geographic structur- 
ing of this karyological diversity, molecular studies re- 
vealed a well-structured organization of genetic diversity of 
mtDNA into two major clades in Atlantic forests (Nogueira 
and Fagundes, 2008). These authors proposed that the di- 
versification of populations in this species was shaped by a 
putative geographic barrier instead of a reproductive barrier 
created by errors in chromosomal pairing. This suggestion 
was apparently supported by the fact that the 2n = 14, 
2n = 15 and 2n = 16 karyotypes were equally distributed in 
both lineages. 

By analyzing the overall geographic distribution of A. 
cursor, Nogueira and Fagundes (2008) proposed that the 
observed genetic variation followed the central-marginal 
model in which the highest genetic variability occurs in the 
center and declines towards the edges of the range. In their 
study, peripheral populations from northern and southern 
lineages showed the lowest intra-populational genetic di- 
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versity whereas central populations from Espirito Santo 
had the highest such diversity. Based on this finding, No- 
gueira and Fagundes (2008) suggested the existence of a 
contact zone between northern and southern lineages some- 
where in the area of the Rio Doce basin (north of Espirito 
Santo), and that new samples from this region should be in- 
cluded in future studies. 

The aim of the present work was therefore to examine 
the partitioning of genetic variability in A. cursor by includ- 
ing mainly samples from Espirito Santo in order to deter- 
mine whether the Rio Doce is as an effective barrier for 
gene flow, thereby stimulating population divergence. For 
this, we included samples from the north and south of the 
Rio Doce basin that had never been tested in previous stud- 
ies. Samples from Sao Paulo (southern lineage) and Bahia 
(northern lineage) were used to represent the two major lin- 
eages previously observed and to help polarize the analysis. 

DNA was extracted from samples of frozen muscle or 
liver obtained from 64 specimens of A. cursor collected at 10 
spawning sites in three states with Atlantic rainforest: Bahia 
(n = 10), Espirito Santo (n = 48) and Sao Paulo (n = 5) 
(Table 1, Figure 1), as described by Bruford et al. (1992). 
Upon extraction, a 1 140 bp fragment of mitochondrial DNA 
containing the cytochrome-b (Cyt-b) gene was amplified us- 
ing the primers described by Smith and Patton (1993). The 
amplifications were done in a 25 uL reaction volume con- 
taining IX PCR buffer (Invitrogen), 1.0mMMgCl 2 ,0.8mM 
dNTPs, 0.2 mM of each primer, 3.0 U of Platinum Tag poly- 
merase (Invitrogen) and 5-50 ng of template DNA. The ther- 
mal cycling conditions were: 94 °C for 3 min followed by 
40 cycles of 94 °C for 30 s, 48 °C for 40 s and 72 °C for 90 s, 
with a final extension step of 72 °C for 5 min. 

The amplicons were purified with either PureLink 
PCR purification kits (Invitrogen) or ExoSAP-IT® (USB 
Corporation), according to the manufacturers instructions, 
and then quantified in a NanoDrop spectrophotometer prior 
to cycle-sequencing. Sequencing used forward and reverse 
primers and was done directly in a capillary ABI3700 ge- 
netic analyser with the Big Dye Terminator protocol (Ap- 
plied Biosystems). Electropherograms were inspected 
visually using BioEdit v. 7.09 (Hall, 1999) and alignments 
were done using Clustal X in MEGA v. 4 (Tamura et al., 
2007). 

ARLEQUIN v. 3 . 1 1 (Excoffier et al. ,2005) was used 
to explore the genetic characteristics and partitioning of nu- 
cleotide diversity. The number of polymorphic sites and the 
haplotype (h) and nucleotide (n) diversity indices were 
computed (Nei, 1987). The appropriate model of nucleo- 
tide sequence evolution was determined using MODEL- 
TEST v. 3.7 (Posada and Crandall, 1998) and AIC parame- 
ters were chosen. Hierarchical analysis of molecular 
variance (AMOVA; Excoffier et al., 1992) done with 
ARLEQUIN v. 3. 1 1 was used to investigate the patterns of 
historical population structure, with the fixation indices 
serving to estimate the proportion of variation within popu- 
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Figure 1 - Map of localities sampled for A. cursor in eastern Brazil and median-joining network of haplotypes. Circle sizes are proportional to haplotype 
frequency and each circle is colour coded according to the haplotype's presence at the sampling locations. Unobserved haplotypes are represented by 
small red circles. Branches are proportional to the number of mutations, except for those with numbers, which represent highly divergent haplotypes with 
more than 10 mutations. Acronyms, geographic coordinates and altitude of localities are: Bahia state, municipality of Una (UNA, 15° 16' S, 39°04' W, 66 
m), Espirito Santo state, municipalities of Aguia Branca (ABR, 18°52' S, 40°49' W, 373 m), Governador Lindenberg (GOV, 19° 16' S, 40°28' W, 220 m), 
Santa Teresa (STE, 19°55' S, 40°36' W, 950 m), Cariacica (CAR, 20°16' S, 40°25' W, 780 m), Viana (VIA, 20°24' S, 40°29' W, 25 m), Domingos Mar- 
tins (DOM, 20°22' S, 40°40' W, 701 m), Castelo (CAS, 20°36' S, 41°1 1' W, 262 m) and Ibitirama (IBI, 20°29'S, 41°43'W, 935-1063 m), and Sao Paulo 
state, municipality of Iguape (IGU, 24°42' S, 47°34' W, 144 m). 



lations (Ost), among the groups (Oct), and among the popu- 
lations within groups (<E>sc). Pairwise cfST values were also 
generated for population comparisons. The amount of on- 
going gene flow between populations was estimated using 
the <S S t values (Hudson et al., 1992). To depict the hierar- 
chical relationship between haplotypes that reveals the 
mutational steps involved in the transition of one haplotype 
to another, and to assess the association between haplotype 
and geography, a haplotype network was generated using 
the median-joining (MJ) algorithm (Bandelt et al., 1999) as 
implemented in NETWORK 4.2.0.1 (http://www.flux- 
us-engineering.com/sharenet.htm). Phylogenetic analyses 
were done using PAUP* 4.0b 1 (Swofford, 2003) with max- 
imum-parsimony (MP). Support for individual clades was 
evaluated using bootstrap re-sampling (Felsenstein, 1985) 
with 1000 replications of random addition and tree bisec- 
tion reconnection (TBR) branch swapping. 

Sixty haplotypes from 64 individuals were detected 
across all populations (Table 1). Of 1140 sites, 205 were 
polymorphic, with high haplotype and low nucleotide di- 
versity indices in all populations. There was high haplotype 
diversity (h) and low nucleotide diversity, such as is fre- 
quently encountered in populations undergoing a popula- 
tion bottleneck followed by demographic expansion, with 
the accumulation of mutations (Grant and Bowen, 1998). 
This feature could explain the high number of polymorphic 
sites and is not uncommon in Atlantic Forest species. 



The MP tree topology (data not shown) and me- 
dian-joining network (Figure 1) showed two major diver- 
gent lineages. The ten haplotypes of the 10 individuals from 
the northern lineage occurred in southern Bahia, whereas 
the 50 haplotypes of the 54 individuals of the southern lin- 
eage all occurred at localities in Espirito Santo (north and 
south of the Rio Doce) and Sao Paulo (Table 1), and were 
separated by at least 30 mutational steps. Pairwise distance 
differences showed no gene flow between northern (BA) 
and southern (ES and SP) populations, and was confirmed 
by the high <Tj S t values between BA/ES and SP populations 
(Table 2; <D ST > 0.63790, p < 0.05). The phylogenetic tree 
topology and haplotype network agreed with the clustering 
of the northern and southern lineages observed in the maxi- 
mum parsimony analysis presented by Nogueira and Fa- 
gundes (2008). 

Akodon cursor had a significant population genetic 
structure of <D S t = 0.51, and AMOVA, which was used to 
assess whether the Rio Doce influenced this group division 
(Table 2), showed that this river did not affect formation of 
the northern and southern lineages. 

The pattern of phylogeographical groups observed in 
A. cursor could be limited by geographical landmarks, such 
as the Rio Doce basin, as observed in other mammals 
(Mustrangi and Patton, 1997; Costa et al, 2000), geckos 
(Pellegrino et al, 2005) and birds (Cabanne et al., 2007), 
although the splitting times among clades do not suggest a 
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common temporal origin for these patterns (Cabarme et al., 
2007). However, since the northern populations in Espirito 
Santo (Aguia Branca and Governador Lindenberg) are lo- 
cated north of the Rio Doce, the network data supported the 
AMOVA results, i.e., that the Rio Doce does not play an 
important role in the divergence of A. cursor lineages and 
therefore that this geographical barrier should not be con- 
sidered as the main cause of the breakage in the dispersion 
of these rodents, contrary to the proposal by Nogueira and 
Fagundes (2008). 

Additional network analysis showed that all five ha- 
plotypes from Iguape (Sao Paulo) grouped as a monophy- 
letic clade (Figure 1). Despite the generally low divergence 
between southern (ES and SP) populations, 4>st analyses 
involving the Iguape (SP) and Aguia Branca or Governador 
Lindenberg (ES) populations showed high divergence (Ta- 
ble 2). 

Although this was no evidence of previous fragmen- 
tation in populations from Sao Paulo, this apparent "mono- 
phyly" may reflect either the lack of sampling of intermedi- 
ate populations (including Rio de Janeiro state, for 
example) or the fact that the SP and ES populations are geo- 
graphically far (-1000 km) from each other in the southern 
lineage. These observations suggest that isolation by geo- 
graphical distance and the history of vicariance (geograph- 
ical location) could play an important role in shaping the 



population genetic structure of the southern lineage, includ- 
ing populations from Sao Paulo and northern Espirito 
Santo. 

It is well-known that sampling the extremes of a pop- 
ulation distributed as cline can yield spurious evidence of 
vicariance. The data for Aguia Branca and Governador 
Lindenberg provide further evidence of this potential isola- 
tion by distance and support the conclusion that in the geo- 
graphically complex Atlantic Forest it is unlikely that forest 
refugia or barrier hypotheses alone will account for general 
patterns of lineage diversification (Thome et al, 2010). 
However, testing the isolation by distance mechanism re- 
quires further investigation. 

Specimens with karyotypes 2n = 14 and 15 were dis- 
tributed across all haplotype network branches (Table 1 and 
Figure 1), thus indicating that karyological data are not as- 
sociated with lineage diversification, a conclusion that cor- 
roborate previous studies in A. cursor (Nogueira and Fa- 
gundes, 2008). 

In conclusion, our data do not support a primary influ- 
ence of the Rio Doce as a putative barrier in shaping the di- 
vergence between the two main mitochondrial lineages of 
A. cursor in the southeastern Atlantic forest. Rather, the 
data raise questions about the influence of isolation by geo- 
graphical distance in modeling the phylogeographical 
structure of the southern lineage of A. cursor. However, 



Table 2 - Genetic structure of the A. cursor populations. Upper half of table: pairwise <1>st (genetic differentiation between populations). Lower half of ta- 
ble: Analysis of molecular variance (AMOVA) to test two hypothetical scenarios for the partitioning of genetic variation using the fixation indices among 
groups (Oct), among populations within groups (Osc) and within populations (Ost). Significant p-values (p < 0.05) are highlighted in bold. UNA (Bahia) 
was used to polarize the analysis and represent the northern lineage in A cursor (see Figure 1). 
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there is uncertainty about effectiveness of the processes in- 
volved, and sampling on a wider scale is necessary to assess 
the importance of distinct isolation mechanisms. 
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